Scaling of the Atmosphere of Self- Avoiding Walks 



A L Owczarek 1 and T Prellberg 2 

' 1 Department of Mathematics and Statistics, 

The University of Melbourne, Victoria 3010, Australia. 
' a.owczarek@ms.unimelb.edu.au 

2 School of Mathematical Sciences 
Queen Mary, University of London 
; Mile End Road, London El 4NS, UK 

t . prellberg@qmul .ac.uk 

*0 ■ June 6, 2008 

a 

^ ■ Abstract 

The number of free sites next to the end of a self-avoiding walk is known as the 
atmosphere. The average atmosphere can be related to the number of configurations. 
Here we study the distribution of atmospheres as a function of length and how the 
number of walks of fixed atmosphere scale. Certain bounds on these numbers can 
be proved. We use Monte Carlo estimates to verify our conjectures. Of particular 
interest are walks that have zero atmosphere, which are known as trapped. We 
demonstrate that these walks scale in the same way as the full set of self-avoiding 
■ walks, barring an overall constant factor. 
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Consider an n-step self- avoiding walk u> = (co , u>i, . . . 
d > 2, and steps having unit length, i.e. \u> i+ i — Ui\ = 1. 



with n + 1 sites u)i G Z d for 



Figure 1: A self-avoiding walk on the square lattice Z 2 with n = 20 steps and an atmo- 
sphere a = 2. 

The number of edges that can be appended to the last visited vertex u n to create an 
(n + l)-step is called the atmosphere of the walk u. Clearly the smallest value of the 
atmosphere is zero, in which case the walk is called trapped. A zero-step self-avoiding 
walk has atmosphere 2<i, and for n > 1 any n-step self-avoiding walk has atmosphere of 
at most 2d — 1. 

We partition the set of n-step self-avoiding walks by the value of their atmosphere. 
Denote by c n the number of n-step self-avoiding walks starting at oo = 0, and by the 
number of n-step self-avoiding walks starting at u = with atmosphere a. 

The subject of this paper is the fraction of n-step self-avoiding walks with fixed atmo- 
sphere, 

(a) 

Pi"' = — , (1) 



and its limiting behaviour as n — ► oo. 

The atmosphere was introduced by Rechnitzer and Janse van Rensburg pQ for self- 
avoiding walks. Since it is well established theoretically, though not proven, that 



c n ~ Aifi n n 



(2) 



they immediately pointed out that the mean atmosphere could be used to estimate the 
connective constant /x and exponent 7 since 

E(a) 
a ac n _ c n+l 



(a) 



(3) 
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Subsequently they have generalised this elegant idea to self- avoiding walks near walls [2] 
and recently to self-avoiding polygons [3J. 

Motivating this paper is recent interest in various subsets of walks that do not trap 
|H E] including prudent walks [6j Ej. The question that naturally arises for each class 
of walks is "Does the numbers of walks in the class scale in the same way as full set of 
self-avoiding walks"; that is, as in equation (J2]) with the same value of \x and the same 
value of 7. 

It is then natural to consider the limits 

(a) 

p(«) = lim — . (4) 
In this paper we demonstate results in the literature 0, [9] can be used to prove 

liminf — > , (5) 

n^oo C n 

and so that if one assumes the limit exists then one has immediately that 

p {a) > . (6) 
Going further and assuming the scaling in equation §2§ we can therefore predict 

c w ~ 4V* 7 - 1 , (7) 

where 

A[ a) = p^Ai . (8) 

Hence walks that trap, that is, those with atmosphere zero will have the same connective 
constant and same value of exponent 7 as all self-avoiding walks. Conversely the same is 
true for walks that do not trap. 

Furthermore we use the Monte Carlo Algorithm flat PERM |10] to provided estimates 
for pffi up to length n = 512 and use the corrections-to-scaling to extrapolate high- 
precision estimates of p( a \ 

Theory Intuitively one would expect that the probability that the end of a long self- 
avoiding walk uj ends at the corner of a cube Q of fixed size, but otherwise does not 
intersect with the cube, is non-zero. Therefore one can append to this walk a short walk 
contained within the cube with fixed atmosphere. Using a related argument, we shall 
prove below that the limes inferior of pn is bounded away from zero. 
Our main result is the following theorem. 
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Theorem 1. Let < a < 2d. Then 

lim inf p^ > , 

n—>oo 

and 

lim sup pffl < 1 . 

n— >oo 

This theorem is very similar to a result concerning tail patterns of self-avoiding walks. 

Definition 2. A self-avoiding walk P = (po,Pi, ■ ■ ■ ,Pk) is called a tail pattern if there is a 
self-avoiding walk u = (uq, u>i, . . . , u n ) such that P = (u n -k, Wn-k+h • • • > w n ). P is called 
a proper tail pattern if for all sufficiently large n there is a self-avoiding walk having P 
as a tail pattern. 

Let c n [P] be the number of ra-step self-avoiding walks with tail pattern P. We have 
the following result, which is taken from [8]. 

Proposition 3. If P is a proper tail pattern then 

l im in f > . 

n— >oo C n 

When applying this proposition for the proof of Theorem [U we find that the main 
difficulty is that the atmosphere of a proper tail pattern P and the atmosphere of a 
self- avoiding walk having P as a tail pattern can be different. 

Proof of Theorem [3 It is sufficient to show that for any atmosphere a with < a < 2d 
there exist a proper tail pattern P a with atmosphere a and the additional property that 
any sufficiently long self-avoiding walk having P a as a tail pattern has also atmosphere 
a. The first inequality then follows as an immediate consequence of Proposition [3J as 
Cn^ > c n [P a ], and therefore 

Cn ) C \P I 

lim inf p£' = lim inf > lim inf — — — > . 

n^oo n— >oo C n n—*oo C n 

The upper bound follows now from the observation that Y^f=o Pn — 1> as 

lim sup p^ = lim sup | 1 — p^ ^ J 

n^oo n— >oo \ , , / 

\ a'f=a / 

<1 + ^ limSU P (-Pn ^) 



-1 - V lim inf < 1 . 



n— >oo 
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The pattern P a is constructed as follows. Consider a self-avoiding walk u which starts at 
the origin, visits precisely 2d — a neighbouring vertices, and ends on the corner of a cube Q 
containing it. Continue wbya self-avoiding walk u' on the set of vertices dQ which have 
distance 1 from Q, and which is Hamiltonian on dQ. The pattern P a which is obtained 
by reversal of steps is a proper tail pattern. It has atmosphere a, and by construction 
the atmosphere of any sufficiently long self-avoiding walk having P as a tail pattern must 
also be a. □ 

Note that the idea of this proof generalises easily to different notions of atmosphere, 
such as a fc-step atmosphere given by the number of (n + fc)-step self-avoiding walks which 
can be grown from an n-step self-avoiding walk. 

Simulation results Backed up by the rigorous results, we now consider results of simu- 
lations for self- avoiding walks on the square lattice Z 2 . Using flatPERM, a flat-histogram 
kinetic growth algorithm, we have grown 10 9 walks at length 512. As a consequence of 
using a flat-histogram method with respect to the atmosphere, we could boost the occur- 
rence of zero-atmosphere walks roughly by a factor of three, with minimal computational 
overhead. 

From Figure [5] it is clear that the quantities p^ approach an asymptotic value quickly 
as n — > oo with only small corrections to scaling. 

Figure [3] shows that these corrections are asymptotically linear in 1/n. We therefore 
conjecture on the basis of our simulations that the limit 

» 



>< a > = lim C " 



n— too c n 

indeed exists. A linear fit can be used to obtain estimates for p^ , and we find that 

p(°) =0.009096(4) (9) 

pW =0.05476(1) (10) 

p (2) =0.22500(2) (11) 

p (3) =0.71114(3) (12) 

for self-avoiding walks on the square lattice Z 2 . 
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Figure 2: A plot of the probability p$ of an n-step self- avoiding walk on the square 
lattice having atmosphere a. Shown are for a = 0, 1, 2, 3 (from bottom to top) on a 
logarithmic scale versus 1/n. We note that there are only small corrections to scaling. 



by the authors. 



References 

[1] A. Rechnitzer and E. J. Janse van Rensburg, J. Phys. A: Math. Gen. 35 L60512, 
2002 

[2] A. Rechnitzer and E. J. Janse van Rensburg, J. Phys. A: Math. Gen. 37 68756898, 
2004 

[3] E. J. Janse van Rensburg and A. Rechnitzer , J. Phys. A: Math. Theor. 41 105002, 
2008 

[4] A. J. Guttmann, J. Dethridge and I. Jensen, Prudent self-avoiding walks and poly- 
gons. In Random Polymers, EURANDOM (June, 2007) 



6 



0.714 



0.713 



0.712 



0.711 



0.05475 



0.05450 



0.05425 



0.05400 




0.001 



0.002 0.003 
1/n 



0.004 




0.2250 



0.2245 



0.2240 



0.2235 



0.00910 



0.00905 



0.00900 



0.00895 




0.001 0.002 0.003 0.004 
1/n 



0.001 0.002 0.003 0.004 
1/n 




0.001 0.002 0.003 0.004 

1/n 
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